Realization of Haldane's Exclusion Statistics in a Model of Electron-Phonon 

Interactions 
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We discuss an integrable model describing one-dimensional electrons interacting with two- 
dimensional anharmonic phonons. In the low temperature limit it is possible to decouple phonons 
and consider one-dimensional excitations separately. They have a trivial two-body scattering matrix 
and obey fractional statistics. As far as we know the original model presents the first example of 
a model with local bare interactions generating purely statistical interactions between renormalized 
particles. As a by-product we obtain non-trivial thermodynamic equations for the interacting system 
of two-dimensional phonons. 



The concept of fractional exclusion statistics was intro- 
duced by Haldane in 1991E1 and since then there has feo^n 
a continuous interest in this subject (see, for exampleoEl) . 
In the ideal case of a 'quantum gas of non-interacting 
particles satisfying fractional exclusion statistics' (the 
Abelian one) the number of states Q a of particles of 
species a linearly depends on numbers of occupied states 
of particles of species b, N b with the same energy: 

dQ a /dN b = -G ab 

where G ab is a matrix of statistical interactions. As has 
been demonstrated by many authors ,□ the thermody- 
namics of such an ideal gas of anyons is described by 
the set of equations 



e a (£;)+G Qfe Tln[l + e- 



E 



(1) 



f/l = -tYM 1 

a,E 



It has been emphasised by different authors, most re- 
cently by Frahmd and Bouwkncgt and Schoutcnsa that 
integrable systems in (l+l)-dimensions represent natu- 
ral realizations of all sorts of exotic Abelian and even 
non-Abelian exclusion statistics (the latter ones are not 
discussed in this paper). This follows from the fact 
that creation and annihilation operators of renormal- 
ized particles in these systems satisfy commutation re- 
lations which explicitly depend on the two-particle S- 
matrix (the Zamolodchikov-Faddeev algebra) being thus 
neither bosonic nor fermionic: 



z+{e x )z+{e 2 ) 

Z a (9 1 )Z b (9 2 ) 
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s a a i(9 1 -e 2 )z+(9 2 )z±(9 l ) 



s a i( 

a,b 

ski 



9 2 )Z b (9 2 )Z a (9 1 ) 
9 1 )Z^(9 2 )Z a (9 1 )+5(9 1 
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Here the operators Z£(0) and Z a {9) create (annihilate) 
the state with rapidity 9 and isotopic number a (recall 
that rapidity is a relativistic parametrization of energy 



and momentum E = m cosh (9, P = msinh#; for simplic- 
ity we shall discuss only Lorentz invariant models, but 
the results have more general validity) . 

Unfortunately, most known integrable models are not 
good pedagogical examples of the effects of exclusion 
statistics because they describe interacting anyons , 
where the interaction derives from a momentum depen- 
dent ifj-matrix. This is why the free energy in these 
theories is described by a complex system of nonlinear 
integral equations which in general, cannot be solved in 
closed form. The ^-matrix in this models becomes purely 
clastic only at high momentum: 



S b b i(0) -> 5 a a 5 b b exp[-iTT S ign9G ab 



(|0|-oo). (3) 



A model where this helds at all rapidities, would be an 
ideal gas of anyons. Examples of such models known 
so far involve long-range interactions (for instance, the 
Haldane-Shastry spin-chainH; see also the discussion inH) . 

In this paper we provide an example of an integrable 
model with local interactions where the two-body S- 
matrix of physical particles has the form (0) from the 
very start. We study a limiting case of an integrable 
model first suggested by FateevQ . The model describes 
a one -dimensional sea of interacting fermions, coupled 
to a two dimensional bath of anharmonic phonons. The 
Hamiltonian density has the following form: 



where 



7i — Tie + T~L P + tii, 



"He 



(4) 



-iR + d x R + \L+d x L + (2 5 /tt) PjRL^L - fi{R+R + L + L) 

describes a one dimensional electron gas of right (R), 
and left-moving (L) fermions ("Thirring model") with 
the chemical potential fj,, and 



1^, 



N-l 



ri— 1 n—1 

+ l( m / /3 )2[ e Wi +e - 2 ^], 



(5) 
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where [P„(x), 4>m(y)} = -iS nm S(x - y), describes N par- 
allel chains coupled to each other by an anharmonic 
nearest- neighbour interaction , ("Toda lattice"), and 



Hi = m(R ] L+ L^R) 



describes the coupling between the two-dimensional 
phonon gas and the Fermi sea. The pure fermionic 
Hamiltonian corresponds to the massive Thirring model 
which is one of the first integrable models known and has 
numerous applications in condensed matter physics. 6 

In the limit (3 << 1 one can expand in (3. In the limit 
N — > oo the phonon spectrum becomes continuous two- 
dimensional with the dispersion law given by 



w 2 (p) = pI + m 2 sin 2 (7rp y /2) 



(6) 



As we shall see this form of dispersion remains robust 
at any /?. We shall be interested in the case when 
tp = fj, — M > when there are charged particles in 
the ground state. It turns out that in this case phonons 
decouple and one can study the remaining edge excita- 
tions separately. In terms of the original fermions the 
corresponding effective action describes one-dimensional 
fermions interacting by a phonon-mediated interaction 
non-local both in time and space. For instance, in the 
leading order in (3 we have 



<< e Wi("<-,x) e Wi(0,0) >> = 



il 2 



47t(t 2 



2)3/2 



(7) 



The exact solution gives the most illuminating results 
for thermodynamics demonstrating that true edge excita- 
tions are ideal anyons. Namely, there are two oppositely 
charged particles realizing Haldane's exclusion statistics. 
Particles with the same charge scatter with phase shift g, 
whilst particles with opposite charge scatter with phase 
shift 7r — g, such that 



G = 



g/n 1 — g/ir 

1 - ah ah 



(8) 



We would like to mention here an amusing feature of 
model (H|): under the transformation 



g -> 7T - g, ((3 -» 4tt//?) 



(9) 



the fermionic part of model (Q) transforms into a non- 
linear sigma model of a single complex bosonic field and 
the phonon part is transformed into a Toda chain with 
N—1 sites (one site shorter) so that the corresponding 
Lagrangian density becomes, according to FateevQ, 

= 1 d t x*d t x - d x \*d x x m 2 2 ^ 
2 1 + (/3/2) VX 2 1X1 6 



N-l 
n=l 



N-2 



{(dt^f - (d x c/> n ) 2 } - {m/(3) 2 c ' 



Therefore a weak coupling limit of the fermionic model 
is dual to the strong coupling limit of the bosonic model 
@ with N—1 sites and vice versa. This duality is re- 
produced in our results serving as a check for correctness 
of the exact solution. The duality was obtained by Fa- 
teev in essentially nonperturbative way and it remains 
unclear at the moment how to relate the fields of model 
© to the original fermions and bosons. n 

Let us recall the relevant results from (0). The spec- 
trum consists of a charged particle (with the correspond- 
ing anti-particle) of mass M and N — 1 neutral particles 
with masses 



M a = 2Msin(7ra/27V),a = 1, 2, ...N - 1. 



(11) 



where M is related to the bare mass m. The latter ex- 
citations are particle-hole bound states. It is very easy 
to establish a connection between these particles and the 
phonons of the Toda lattice: the above spectrum is re- 
produced in the limit (3 — > when the phonons become 
harmonic (see Eq. (|^) . It is a well known fact in the theory 
of Toda models that the exact spectrum coincides with 
the one CjbtaiHed in the .harmonic approximation (see, for 
example,B,El,ll3; see alscO for general discussion of ther- 
modynamics of theories with elastic S- matrices). The 
scattering is purely elastic (all S matrices are diagonal in 
isotopic indices). 

The Thermodynamic Bethe Ansatz (TBA) equations 
for models with-elastic scattering has the following stan- 
dard form (seeO for a review). The free energy of the 
system of length L is given by 

F/L = -TJ2^ /"d0cosh#ln[l + e- £ "W/ T ] (12) 

a 

where the excitation energies are determined from the 
following nonlinear integral equations: 

Tln[l + e e "^/ T ] - (aS nm + G nm ) * 7Tn[l + e'^' T ] 

= m„cosh# (13) 



G n m(9) — ~ 



1 dlnS nm (0) 
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where a = 1 for fermions and a = for bosons (that is 
for the neutral bound states) and the star denotes con- 
volution: 



G * f(6) 



G(6 - e')f(e')de' 



From RefS we obtain the following expressions for the 
Fourier transformations of the kernels: 



Gabiu) 



4sinh(5cj/2A) sinh[(7r - g)u)/2N] 
sinh(7rw/2A) 



(10) 



cosh[(iV — max(a, 6))7rw/27V] sinh[(mm(a, fyiruj/IN] 
cosh(7rw/2) 
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G„ J+ (w) = G„,_(w) 



2 sinh(.gcj/2iV) sinh[(7r - g)u/2N] sinh(a7rcj/2iV) 
sinh(7rw/2iV) cosh(7ra>/2) 



G++M = G — (w) = 



4sinh(c/w/2iV) cosh[(7r(7V - 1) + g)w/2N) 



sinh(7ru;/2A0 



cosh(7ra;/2) 



G+- (w, .9) = G++ (w, 7r - g) 



(14) 



where "+" and "-" stand for particle and anti-particle 
and "a,b" label phonon branches. 

In the limit N — > oo the system ([13]) undergoes a dras- 
tic simplification. In particular, we get 



G++M = -g/ir,G+-(u) = -(l-g/ir) 



(15) 



that is the kernels in the TBA for the ^-particle and 
antiparticle become delta functions. 

Another important fact is that in this limit the 
phonons decouple from the TBA equations for "+,-" and 
can be analyzed separately. To pass to the N — ► oo limit 
we introduce a continuous variable y = a/N and replace 
summation over a, b by integration over y: 



E -> N I d y 

Jo 



Eqs.(p"3|) become 

rin[l + e e ±W/ T ] - (1 - gMTln[l + e~ e ±W/ T ] 

+ (1 -g/^)T\n[l + e-^' T ] 

= M cosh ± /u 

+ T d V J d0 ' K + iv, - M 1 + e~ £( ^' )/T ] (16) 

where the Fourier transform of the kernel is 

K+(y,u>) = 2g(l — g/ir)uismh(yTru)/2)/ cosh(7ru;/2) 
and 

Tln[l + e<^)/ T ] 
-T J dy 1 J d9'K(y, y'; 6 - 9') ln[l + e -<v'#')/ T ] 
= 2M sin(7ry/2) cosh 6 (17) 
where 

K(y,y';uj) = 4g(l -g/-K)ux 



cosh[7r(l — max(y, y'))u>/2] smh[irmin(y, y')u)/2] 
cosh(7rw/2) 

Note that terms with e± are absent in the last equation, 
as expected in the large- N limit where the fermions pro- 
vide a small boundary effect on the phonons giving a 
0(l/7V)-correction to their thermodynamics. 

Eq.(jl^) is further simplified in the most interesting 
limit 



M/N < T < M 



(18) 



where only the region y <C 1 contributes to the thermo- 
dynamics. Defining the new scaling variable £ = nyM/T 
we get 



hi[l+e e ^/ T ]+ 7 



de 2 



dC'mz7i(^,^)ln[l+e" e(s/ '' f,)/T ] 



£ cosh 9 



(19) 



where 7 = 2g(n - g){T/irM) 2 . 

This equation together with the expression for the free 
energy 



Fphon/NL 



rp3 p OO 

/ d^d6»cosh6»ln[l + e- e( «' e)/T ] (20) 
\tt 2 M J 



describes the continuous limit of the (2 + l)-dimensional 
Toda array. As far as we know such limit has never been 
described before. It follows from Eqs.(|n],|2(]) that tem- 
perature enters into e/T only through 7 and therefore 



T 3 



Fphon/NL = -— Y [2 9 {-k - g)(T/nMf 



(21) 



where Y(-j) is a smooth positive function allowing reg- 
ular expansion in 7. This result is valid in all powers 
of the bare coupling constant (3; being dependent on 
g(n — g) = 47r 3 (47r//3 + 0)~ 2 it is manifestly self-dual 
under transformation (^). 

Eqs.( |l9| , p0|) can be studied by iterations in 7. For this 
purpose it is convenient to introduce the quantity 



ln[l 



such that these equations become 

-)2 



d 2 f°° 1 
E ^ e )-^QffI J o d^min(e,Ola[l-e- B « = £cosh0 



Fphon/NL 



2tt 2 M 



d££d6>cosh6>ln[l -e 



-E(£,6)i 



(22) 



In this formulation it is absolutely clear that in the limit 
g — > 0, 7r the Toda array becomes equivalent to the free 
(2+l)-dimensional bosonic field and that the phonon free 
energy allows a regular expansion in 7 ~ T 2 . Indeed, the 
substitution p x = M£smh6,p y — M£ brings Eq.(p2|) to 



3 



the familiar form of the free energy of two-dimensional 
acoustic phonons. The first terms in the expansion in 7 
are 



1 



Y(j) = -C(3) + 



2tt 2 



ydy 



1 



x dx 
e x - 1 



0( 7 2 ) (23) 



Further we shall be interested in the case when chem- 
ical potential overcomes the gap M and there is a finite 
density of charge in the ground state. In this case the 
fermionic free energy is expanded in powers of (T/ep) 
where ep ~ fi — M is the Fermi energy. In comparison 
with fcrmions two-dimensional phonons give extra small 
factor T/M. Different temperature dependences of the 
fcrmion and phonon sectors allow one to consider them 
separately. 

To illustrate this point let us consider T the last term 
in the right hand side of Eq. . At T <C M it becomes 



and can be treated as a thermal renormalization of the 
particles dispersion. We shall omit it together with renor- 
malization of the fermion free energy coming from polar- 
ization of the host - the latter term being of order of T 3 
in comparison with T -contribution coming from the ID 
particles. 

After these simplifications the effective TBA for the 
fermions become 



-(l- S /7r)Tln[l + 



( 3 /7r)Tln[l + e- e + (e)/T ] 



McoshO-n (25) 



e_(0) + (. 9 /7r)Tln[l + e- £ -W/ T ] 
+ (l-.g/7r)Tln[l + e- e+(9)/T ] =Mcosh0 + ,u (26) 

with the free energy equal to 
F/L = -T^2 — [ d0cosh01n[l + e- e «W/ T ] (27) 

a=± n J 

These equations coincide with the general equations for 
particles with exclusion statistics (|l|) for the particular 
choice of the matrix G a b given by (||) (the fact that exclu- 
sion statistics exists in this model has been first pointed 
out inE3). As an illustration we give their explicit solu- 
tion for g/n =1/2 and T << M: 



of particles obeying exclusion statistics. As a byproduct 
we have written down equations for the free energy of 
two-dimensional model of anharmonic phonons. As we 
have pointed out in the introduction, one purpose of this 
excersise is pedagogical—It has been pointed out many 
times (for instance, inaji3) , that the ideal anyon gas in 
one dimension is just a generalization of Luttinger liquid 
for non-linear quasi-particle dispersion thus having many 
properties common with Luttinger liquids. On the other 
hand, analysis of such physical processes as transport 
looking rather complicated in the standard bosonization 
appmach is greatly simplified in the anyon gas formalism 
(seeE 3 ]). We believe that the model discussed in this paper 
will be valuable for such analysis. We also emphasise that 
the model discussed is the first example of a solvable 
model with electron-phonon interaction. 
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F/L = —T— / dtfcoshflln 

7T 



1 



-eo(S)/T 



+ jl e -2e (e)/T + 1 



where e (8) — M cosh. 8 — fi. 

In summary, we have found a relatively simple and 
realistic model which turns out to be a model of ideal gas 



(28) 
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